Macroscopic Many-Qubit Interactions in Superconducting Flux Qubits 
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Superconducting flux qubits are considered to investigate macroscopic many-qubit interactions. 
Many-qubit states based on current states can be manipulated through the current-phase relation in 
each superconducting loop. For flux qubit systems comprised of N qubit loops, a general expression 
of low energy Hamiltonian is presented in terms of low energy levels of qubits and macroscopic 
quantum tunnelings between the many-qubit states. Many-qubit interactions classified by Ising 
type- or tunnel-exchange interactions can be observable experimentally. Flux qubit systems can 
provide various artificial-spin systems to study many-body systems that cannot be found naturally. 

PACS numbers: 74.50.+r, 85.25.Cp, 03.67.Lx 



Introduction.— The last decade has seen rapidly de- 
veloping advanced material technologies that make it 
possible to investigate previously inaccessible quantum 
systems for quantum information and computation in 
solid-state systems. Especially, coherent manipulation of 
quantum states in tunable superconducting devices has 
enabled to demonstrate macroscopic qubits P, 0, 0] and 
entangled states of qubits @, H, Experimentally, it 
has been shown that, in terms of pseudo-spins, differ- 
ent types of exchange interactions between two artificial- 
spins such as an Ising interaction for charge qubits [4| 
and flux qubits [!, 0] and an XY interaction for phase 
qubits [||, can be realized and controlled by the system 
parameters. 

It is believed that electrons are interacting in a pair. 
The interaction is called two-body interaction. Nor- 
mally, in a number of spin systems such as spin chains 
and lattices, the two-body interactions in the spin pairs 
reveal rich many-electron physics. Understanding the 
many-electron effects is one of the most important re- 
searches in condensed matter physics. Additionally, in 
a strongly correlated electronic system, a low energy 
spin Hamiltonian can involve more than three spin in- 
teractions [H, 0, Such multiple spin interactions 
are known to play a significant role for quantum phase 
transitions. However, multiple artificial-spin interactions 
are not yet investigated, although artificial-spin exchange 
(two-body) interactions are demonstrated in different 
types of solid-state qubit systems. This work aims to dis- 
cuss, in a general framework, how artificial-multiple spin 
interactions are possible and realizable in superconduct- 
ing qubit systems. In fact, flux qubit systems are shown 
to have an intrinsic property which is multiple artificial- 
spin interactions. Accordingly, flux qubit systems enable 
to study various artificial-spin systems corresponding to 
many-body systems unlikely found naturally. 

In a superconductor, the macroscopic wavefunction 
can be written by ip(r) = y/n* e liptyr \ where n* and ip(r) 



are the density and phase of Cooper pairs, respectively. 
ip(r) describes the behavior of the entire ensemble of 
Cooper pairs in the superconductor. The supercurrent 
density in electromagnetic field is given by 



m* 



ftVp(r) - g*A(r) , 



(1) 



where q* and m* are respectively the charge and mass 
of Cooper pairs. Then, the current states of flux qubit 
loops are influenced by the variations of the phase y(r) 
across Josephson junctions and the vector potential A(r). 
A change of current state in a qubit loop results in a 
change of current states in other qubit loops because (i) 
the change of Josephson junction phases in superconduct- 
ing loops coupling qubits mediates the change of the cur- 
rents states of all qubits and (ii) the circulating current in 
the qubit produces the induced magnetic flux that influ- 
ences on all other qubits. In experiments, several ways 
to make two- or four-flux qubits interacting have been 
employed. Disconnected supconducting loops, as the in- 
direct way, are coupled inductively 0, 0] by means of 
the induced magnetic flux. Other direct ways are to in- 
troduce connecting superconducting loops [111 [H, EH, 
which is called phase coupling. Consequently, many flux 
qubits defined by current states can interact all together, 
which can be observable in experiments. 

We present a general expression of A-qubit Hamilto- 
nian describing low energy physics. The Hamiltonian is 
determined by the low level energies and the tunneling 
amplitudes between A-qubit states in the flux qubit sys- 
tems. We define two types of many-qubit exchange inter- 
actions originating from the energy differences of many- 
qubit states and the macroscopic quantum tunneling be- 
tween the states. Further, it is shown that a specific 
coupling scheme enables to map flux qubit systems into 
many-body systems. 

Model— We consider a general model including the in- 
ductive and phase coupling ways. The N flux qubit sys- 



2 



terns are composed of N qubit loops with N' loops con- 
necting the qubit loops. Primed (unprimed) indices will 
indicate qubit (connecting) loops. The charging energy 
of Josephson junctions in the N(N') qubit (connecting) 
loops is given by 



N' 



i' — l a.' 



'Pi' a' 



(2) 



where C(C) are the capacitance of the Josephson junc- 
tions in the qubit (connecting) loops. $o = h/2e is the 
unit flux quantum. a(a') rely on the number of Joseph- 
son junctions in the qubit (connecting) loops. The in- 
ductive energy is given by 



N N,N' 



n L = 



HJ=1 



i, i'—l 



N' 1 



(3) 



i',j' = l 



where -£»(-?,'/) are the circulating currents in the qubit 
(connecting) loop — Lg is the self-inductance 

for the qubit loop i. For i ^ j, LW) is the mutual in- 

(i) 

ductance between t he q ubits i and j. L K is the kinetic- 



inductance [111 Il4l |15[ | in the qubit loop i. Similarly, 
L'g , , and L ttyl ■? ) are denoted for the connecting 
loops. •* is the mutual inductance between the qubit 
loop i and the connecting loop i' . Finally, the Josephson 
energy of the junctions is given by 



Hj = 



N 

EE 

= 1 a 



N' 



2^sin 2 ^ + ^^2^;sin 2 ^, (4) 



i' — l a' 



where Ej's are the Josephson energy of junctions in the 
qubit and connecting loops. 

Fluxoid quantization.— By integrating Eq. |T]) along 
the closed path in the i-th loop, the fluxoid quantization 
gives the boundary conditions, 



fii 



(5) 



where iff is the phase across the Josephson junction a, 
m is an integer, and /j = f^ t + f.^ d consists of an exter- 
nal and induced magnetic fields, i.e., / e *i = &i/&o an d 

/ffi = EU L{l3)l ^^ + Eti^ l ' ] n : /^o- Similarly, 
the boundary conditions in the connecting loops can be 
given. From the boundary conditions, the total energy 
can be reexpressed as a function of the phases, {ipi} and 
their time derivatives, 

N-qubit Hamiltonian.— The number of Cooper pairs, 
n, and the phase of wavefunction, ip, are non-commuting 
variables, i.e., [ip, n] = i, such that the canonical mo- 
mentum, P v , can be introduced as P v — nh — —ihd v 



[lq |. where n — q/2e with the charge from the Joseph- 
son relation, q = C(Q$/2ir)(p. When the charging energy 
is much smaller than the Josephson energy, the phase 
is well-defined while the number is strongly fluctuating. 
The charging energy TLc({<Pi\) plays a role of kinetic en- 
ergy for a particle in an effective potential defined by 
U({ip i }) = H L ({<p i }) + Hj(W i }). 

In the three-Josephson junction qubit loops (a G 
{a,b,c}) with Ej? — Eji and ip\' c — ipi, the effective 
potential has the 2 N local minima corresponding to the 
2 N basis, {|mi, • • • , mpf)} , of the N qubits with m; =| 
and J, for i = 1, • • • , N. The values of {(fii} at the local 
minimum corresponding to the state |mi, • • • , rnjy) are 
denoted by W° mi ... mN }. Then, {^° mi ... mjv } determines 
the current state of flux qubit i by the current-phase re- 
lation. 

In the low energy limit, one can employ a tight- 
binding approximation in which the 2^ states of N 
qubits correspond to 2 JV -lattice sites. In the 2 N ba- 
sis, {|toi, • • • , tun)}, the low energy N qubit-Hamiltonian 
matrix can be written as 



H 



N 



Ji- 



E c h . 

,3N£{0,x,y,z} 



' N > 



(G) 



where (T (cr x ' v,z ) are the identity (Pauli) matrices. The 
coefficients are obtained by 



'JV 



H 



N 



(7) 



The diagonal components of the Hamiltonian matrix 
are the level energies, £ mi .... mjv , at the local minima, 
{ l Pi-m 1 ... mN }- The level energies are given by 



fc N 

2 / ^ i,mi,-,m» 
i=l 



+ U(Wl mi ... mN }), (8) 



the 



characteristic oscillating frequencies are 
' U({<pi})\{ v o y with an effec- 



(i) 



where 

, .2 

tive mass A'U — (f^)~ C^g and effective capacitance C^g 
in the harmonic oscillator approximation [17| . 

Generally, the macroscopic tunneling processes be- 
tween any two many-qubit states are possible due to the 
quantum fluctuation originating from the kinetic energy. 
The off-diagonal components are the macroscopic quan- 
tum tunneling amplitudes, i.e., 



t : |, 



' JV / 



mi, 



,m N ) 



(9) 



for the tunneling between the two states, |m^, • • • ,m' N ) 
and |mi, • • • , tojv) . The tunneling amplitudes can be 
calculated by the well-known numerical methods such 
as WKB approximation, instanton method, and Fourier 
grid Hamiltonian method [HI]. The tunneling process, 
I TTT ' ' ' T) |4TT ■ ■ ■ T)> describes the first pseudo- 

spin flip. Such a tunneling process that describes only 
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one pseudo-spin flip among the N qubits is called single 
qubit tunneling, t±. If the N qubits are flipped for tunnel- 
ing, the tunneling processes can be called N-qubit tunnel- 
ing, t N , e.g., |J.tt • • • I) ^=> ITU •••?)• Normally, single 
qubit tunneling amplitudes are much larger than other 
multiple qubit ones. However, when a multiple qubit 
tunneling amplitude is larger than single qubit one, the 
multiple qubit tunneling processes can play an impor- 
tant role in determining the property of eigenstates of 
the system 1|| 20j |. 

Many-qubit interaction.— Actually, Eq. ^ describes 
any N qubit system including all types of many-qubit 
interactions. Let us expand the low energy N qubit- 
Hamiltonian matrix in terms of qubit interactions; 



i<j<k 



(ijk) 
3 



-H 



(l-N) 



N 



(10) 

where Hq = (l/2 JV )Tr [Hn] and the qubits are described 
by H^f 1 = e% <j\ + erf with the energy difference 2ej 

(i) 

and the tunneling amplitude t\ between the two states 
of the qubit i. Qubit interactions are denoted by two- 
qubit interactions i three-qubit interactions H^® , 
and so on. Then, the JV-qubit interaction is presented by 



H 



(l-N) 
N 



c, 



N 



(11) 



We define the N-qubit exchange coupling constant as 



7 W _ r - 



1 



— — , A -Tl [(T 1 



cr* N Hn] , (12) 



which has a form of Ising type exchange interaction for N 
qubits. For other terms of the iV-qubit interaction, the 
coefficients of the terms can be called N-qubit tunnel- 
exchange coupling constants, e.g., 



J : 



(JV) 



G x 



1 



r z N Hn 



since the off-diagonal components of the Hamiltonian ma- 
trix result from the hopping (tunneling) between the sites 
(states). 

Two qubit systems.— For two qubit systems, the two- 
qubit interaction is given by H. 

cr 2 , where J : 



(12) 



(2) 



(i| + t b 2 )/2, 4 2 y } = (q - t»)/2, and 



2^je{x, v , z } J jj °V 



The two-qubit tun- 
describe the tunnel- 



J ( zV = (E n -E n -E u +E u )/4. 
neling amplitudes, (i) t 2 and (ii) t\, 
ing processes, (i) |||) |JJ) in the parallel pseudo-spin 
states and (ii) ||J.) <J=^> |J,|) in the anti-parallel pseudo- 
spin states. As expected, the exchange coupling constant 

(2) 

Jz Z is the energy difference between the parallel and 
anti-parallel pseudo-spin states. The two-qubit tunncl- 
ings contribute to the pseudo-spin exchange interaction. 

( 12) 

Then, H 2 has a form of XYZ model for two-pseudo 
spins. t 2 <C t 2 gives an XXZ pseudo-spin model and, for 



r(2) 
J zz 



0, 



E 



rr ■ 



E 



U 



E+\ +E\-t, an XY pseudo-spin 



model. For t 2 ' b <C J^z , H 2 12 ^ becomes an Ising pseudo- 
spin model. This shows that various types of pseudo-spin 
models can be realized by manipulating the system pa- 
rameters. 

Three qubit systems.— Next, for comparison, let us 
consider a two-qubit interaction of three qubit system 



given by H^ 2) = J2 

J zxQ cr l 
t( 2 ) _ (Ta 
J yyO - \ l 2 



je{x,y,z} J ijO (J l K 
r(2) 



" JxzO (T l ' 



crl where J^ = -(t a 2 + t§)/4 - (tg + *g)/4, 



*§)/*+ GS-t§)/4> and J« = {E m -E m - 
Em+ E m)/% + {Em-E ni -E iu + E m )/&. Here, i(t) 
denote the two-qubit tunnelings for the up(down) state 
of the third pseudo-spin. Compared to the two-qubit in- 
teraction in two-qubit systems, interestingly, there are 

(2) 

the two extra tunnel- exchange coupling terms, J xz0 and 

(2) 

Jzx0> mediated by the single-qubit tunnelings. 

(123) 

In the three-qubit interaction H 3 , the three-qubit 

( 3) 

exchange coupling constant is given by J Z zz — (^TTT — 

E ll1 _ ^ITT + E ll\)l% _ ( E 1U - ^TU - E IU + E lll)/&- 
Also, the single- and two-qubit tunnelings as well as 
the three-qubit tunnelings give rise to the three-qubit 
tunnel-exchange coupling constants. Especially, if the 
three-qubit tunnelings are stronger than the two-qubit 
tunnelings, the ground state can be in a Greenberger- 
Horne-Zeilinger (GHZ) and if the two-qubit tunnelings 
are stronger than the three-qubit tunnelings, a W-state 
can be generated in an excited state [201 ]. 

Multiple qubit systems.— To explore many-qubit in- 
teractions explicitly, let us consider a specific multiplc- 
qubit system in Fig. [T] (a). For simplicity, the induc- 
tances are assumed to be very small and then the in- 
ductive energy can be negligible. The boundary con- 
ditions for (i) the qubit loops and (ii) the connecting 
loop are reduced to (i) 2ip i + ipf = 2ir(ni — f^ t ) and (ii) 



U-- 



= 2im' 



The effective potential is given as 



4Ej sin 2 + 2E a j sin 2 n(rH 



Ai) 
J cxt 



2E' J sm 2 iT ri -Y^ = i{n. 



f (i) 

J cxt 



Pi/*) 



For the four qubit system (JV = 4), we plot the ex- 
change coupling constants as a function of / = / c ^{ 
and E'j in Fig. [T] (b) and (c), respectively. At the 
co-resonance point, / = 0.5, the three-qubit interac- 
tion disappears while the two- and four-qubit interac- 
tion strengths reach their maximum values in Fig. [T] (b) . 
The sign of the three-qubit interaction is changed from 
negative for / < 0.5 to positive for / > 0.5. As E'j 
increases, the two-, three-, and four-qubit interactions 
increase monotonically in Fig. [T] (c) . 

Interestingly, the four-qubit interaction is stronger 
than the two- and three-qubit interactions. That is, 



J: 



(4) 



3 J 



(2) 



Also, the three-qubit interaction can 
be stronger than the two-qubit interaction for a certain 
applied magnetic field. This result seems to be coun- 
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FIG. 1: (Color online) Top: (a) A N flux qubit system with 
one connecting loop. The N superconducting loops are con- 
nected by a connecting loop interrupted by a Josephson junc- 
tion, E'j. In each qubit loop, the diamagnetic (paramagnetic) 
current states assigned by ||) (|T)) 5 are superposed, which 
makes the loop being regarded as a qubit. (oppositely 
®) denote the directions of the applied and induced magnetic 
fields, fi = $i/$o, in the qubit loop i. h{I') stand for the cur- 
rents in the qubit i (connecting) loop. Ej's are the Josephson 
coupling energies of the Josephson junctions in the connect- 
ing and qubit loops. The fluxoid quantization in the connect- 
ing loop gives rise to the boundary condition connecting the 
phases, <pi, across each Josephson junction. Both the mutual 
inductances and the fluxoid quantization make it possible to 
realize many-qubit interactions in the N flux qubit system. 
The many-qubit interactions are defined in the text. Bottom: 



j: 



izooi m t ne four qubits (N = 

(*) 



T W 7(3) 

J zzz0' 



4) as a function of (b) the ap- 
plied magnetic field / = f™. (i = 1, • • • ,4) for E'j = O.lbEj 
and (c) the Josephson energy E'j for / = 0.45. Other param- 
eters are rn = n' = and Ej = Ej(a 6 {a, b, c}). 



tcrintuitive. However, for an TV qubit system, the re- 
sult can be understood from Eq. ([5]) and the boundary 
condition of the connecting loop, without the assump- 
tion, L' K I'/<!> = n'- (1/2tt)(^ + E,=i Vi) ~ fU where 



f! nd = LT/$c 



J2iLi ^-mIi/^o with the mutual induc- 
tance C^j. When one superconducting loop couples all 
qubit loops, all qubits are interconnected through the 
effective flux f e g = (l/2n) J2iLi ft as weu as fLd- Nor- 
mally, the induced flux is much smaller than the effective 

flux, i.e., /i n d'/ind < Jeff, so that much stronger many- 
qubit interaction for the effective flux than for the in- 
duced magnetic flux can be expected. Therefore, if the 
./V-qubit interaction is much stronger than other quit in- 



teractions 
systems 



H 



21 



A' 



H 



(l-JV) 
N 



can map higher dimensional 



We also considered two more models. (i) For N 
qubits inductively coupled without any connecting loop, 



multiple qubit interactions are intrinsically involved but 
their strengths are very weak, for instance of N = 4, 
Jzzzz ~ 10~ 6 </^oo m tn e parameters of Ref. 0- If the 
two-qubit interactions are much stronger than other mul- 
tiple qubit interactions, 7ijy ~ 2»<J•^2*'^• Then, an 
N qubits inductively coupled can be a many-body sys- 
tem in which one artificial-spin interact with all other 
artificial-spins by the two-body interactions, (ii) For the 
model of Ref. [ll( . the multiple-qubit interactions be- 
have similarly with the model of Fig. [T] (a). In the same 
parameter values with Fig. [T] (b), however, this model 
gives Jzzzz ~ 0.17J^q . The two models show that the 
four-qubit interaction is smaller than the two-qubit in- 
teractions for the four qubit systems. In general, hence, 
many-qubit interactions are dependent on specific exper- 
imental setups and varying the system parameters. Var- 
ious types of artificial-spin systems can be prepared in 
flux qubit systems. Therefore, it is possible to explore a 
many-body system realized in flux qubit systems. 

Summary.— We investigated many-qubit interactions 
in superconducting flux qubit systems. There are two 
types of many-qubit exchange interactions. One is sim- 
ilar with the Ising spin interaction, the other types of 
exchange interactions are due to macroscopic quantum 
tunnclings between the many-qubit states. Various types 
of many-qubit interactions can be realized experimentally 
in flux qubit systems. Moreover, an experimental setup 
can be provided to study many-body systems that can 
be mapped into many-flux qubit systems. 

Cho acknowledges the support from the Australian Re- 
search Council. 
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